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ABSTRACT 



Massive stars form deeply embedded in dense molecular gas, which they stir and heat up and ionize. During an early phase, the 
ionization is confined to hypercompact Hn regions, and the stellar radiation is entirely absorbed by dust, giving rise to a hot molecular 
core. To investigate the innermost structure of such high-mass star-forming regions, we observed vibrationally excited HCN (via the 
direct ^-type transition of V2=l, A7=0, 7=13, which lies 1400 K above ground) toward the massive hot molecular cores G10.47+0.03, 
SgrB2-N, and SgrB2-M with the Very Large Array (VLA) at 7 mm, reaching a resolution of about 1000 AU (0.1"). We detect the 
line both in emission and in absorption against Hn regions. The latter allows to derive lower limits on the column densities of hot 
HCN, which are several times 10 19 cnr 2 . We see indication of expansion motions in G10.47+0.03 and detect velocity components 
in SgrB2-M at 50, 60, and 70 km s relative to the Local Standard of Rest. The emission originates in regions of less than 0.1 pc 
diameter around the hypercompact Hn regions G10.47+0.03 Bl and SgrB2-N K2, and reaches brightness temperatures of more than 
200 K. Using the three-dimensional radiative transfer code RADMC-3D, we model the sources as dense dust cores heated by stars in 
the Hn regions, and derive masses of hot (>300 K) molecular gas of more than 100 solar masses (for an HCN fractional abundance of 
10~ 5 ), challenging current simulations of massive star formation. Heating only by the stars in the Hn regions is sufficient to produce 
such large quantities of hot molecular gas, provided that dust is optically thick to its own radiation, leading to high temperatures 
through diffusion of radiation. 

Key words. ISM: molecules - ISM: structure - ISM: clouds - Stars: formation 



1. Introduction 

The formation of massive stars occurs deep inside dense molec- 
ular cloud cores. The high luminosity of massive (proto)stars 
stirs and heats up the surrounding gas, leading to the evapo- 
ration of ice mantles around dust grains and to high degrees 
of excitation of the molecules. This stage, just prior to ioniza- 
tion and destruction of the molecular gas, is characterized by 
strong line and dust emission fr om a compact r egion, and is 
called a hot molecular core (e.g. ICesaronill2005l) . Since these 
sources are short-lived, they are rare and hence far away, which 
hampers the determination of their physical structure (density, 
temperature, velocity, molecular abundances). Nonetheless, such 
knowledge is a key informa tion for understand i ng the process 
of massive star formation (IZinnecker & Yorkel |2007 | ) and for 
the interpretation of molec u lar line surveys (e. g. Schilke et alj 
2006; iBelloche etafll2007t iBerginet al.l l201Qh . and has been 
const rained by radiative transfer modeling of both single-dish 
(e.g. lHatchell & van der iak|l2003t iRolffs et alJl20Tlt) and inter- 
ferometry observations (e.g. lOsorio et al. 2009). 

Especially interesting is the hottest molecular gas, which is 
closest to the heating sources. A good tracer is vibrationally ex- 
cited HCN (V2=l), whose rotationally excited levels lie ~ 1000 
K higher in energy than the ground vibrational state, and which 
has transitions in the cm-wave regime, whose radiation is not 
obscured by dust. These so-called direct {-type transitions are 
intrinsically weak, but have been detect ed in Galactic high-mass 
star-forming regions (iThorwirthlEOQlT) as well as toward the 



proto-planetary nebula C RL61 8 ( Thorwirth et al.l2003bb and the 
starburst galaxy Arp220 (ISalter et al.ll2008h . 



Obviously, this hot gas is very compact, and high spatial 
resolution is needed to investigate it. In this paper, we present 
VLA observations (with among the highest resolution ever ob- 
tained for thermal molecular line emission) of three massive hot 
cores which have strong HCN direct {-type lines (as measured 
with the Effelsberg 100-m telescope, Sect.l3.4l >. The sources are 
G10.47+0.03 (e.g. IWvrowski et al.ll 1999h . which lies at a dis- 
tance of 10.6 kpc (Pandian et al. 2008) and has a luminosity of 
about 7 x 10 5 L Q dCesaroni et alJl2010h . and SgrB2-N and -M 
near the Galactic center at 7.8 kpc distance dReid et al.l 12009). 
which have total luminosities of about 8.4 x 10 5 and 6.3 x 10 6 
L , respectively (Golds mith et all l992). All sources are sites of 
very active high-mass star formation and contain several Hn re- 
gions. 



After describing observations and data reduction (Sect. |2}, 
we present the results (Sect. |3) and the radiative transfer model- 
ing (Sect.|4]i, discuss the implications (Sect. [5]l and draw conclu- 
sions (Sect. [6}. 
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Table 1. Observational summary 



Source 


Pointing center 


Continuum beam" 


Line beam' 3 


rms 1 km s 16 


rms continuum 6 


rms integrated line 6 




RA(J2000) Dec(J2000) 


(";") 


(";") 


(mJy/Beam) 


(K) 


(Kkms- 1 ) 


G10.47+0.03 


18 08 38.236 -19 51 50.3 


0.12x0.06; 82 


0.14x0.1;69 


2.3 


45 


900 


SgrB2-N 


17 47 19.902 -28 22 17.8 


0.12x0.08; 60 


0.15 x0.11;51 


3.0 


75 


1300 


SgrB2-M 


17 47 20.202 -28 23 05.3 


0.12x0.08; 63 


0.15x0.11;51 


3.3 


95 


1400 



Notes. (a) Different weighting schemes were used for continuum and line images. Major and minor axis and position angle are given. (6> rms noise 
in 1 km s channel, continuum and integrated line maps (Figs. [2JH] and[6j. The noise is evaluated over a box from -2 to -1" in both RA and Dec 
relative to the phase center. 



2. Observations and data reduction 

2.1. VLA 

With the National Radio Astronomy Observatory's (NRAO'sfl 
Very Large Array (VLA) in its BnA configuration we observed 
three massive star forming regions in the direct ^-type line of vi- 
brationally excited HCN at 40.7669 GHz (project AR687). This 
transition connects the two levels of the V2=l, 7=13 doublet, 
whose energies lie 141 1 and 1413 K above the ground state, re- 
spectively. 

Each of the sources was observed during LST 15-21; 
G10.47+0.03 on 2009 January 30, SgrB2-M on January 31, 
and SgrB2-N on February 1. Bandpass calibration was done on 
3C279 and 3C454.3, and primary flux calibration on 3C286. 
The observations were performed in fast-switching mode with 
150s on the source and 30s on a phase calibrator (1755-225 for 
G 10.47 +0.03 and SgrA* for SgrB2), yielding an on-source inte- 
gration time of about 4.5 hours per source. For amplitude cal- 
ibration, 1733-130 (G10.47+0.03) or SgrA* (SgrB2) was ob- 
served every 20 minutes, and reference pointing was done every 
hour on these sources. 

For G10.47+0.03, we centered one IF with 6.25 MHz band- 
width at 40.762 GHz (shifted from the rest frequency due to 
the apparent source motion) and one with 25 MHz bandwidth at 
40.78 GHz to increase the continuum sensitivity. For the SgrB2 
sources, the line is broader, so we used two overlapping IFs with 
6.25 MHz bandwidth centered at 40.761 GHz and 40.765 GHz. 
The channel width is 0.72 km s _1 . 

Calibration and editing was done in AIPS: Bad data points 
were flagged; phase, amplitude, and flux were calibrated us- 
ing the continuum data of the calibrator s ources, and band - 
pass calibration was applied. In MIRIAE0 (S ault et all 11991 . 
Doppler correction was done and the channels were averaged to 
a width of 1 km s _1 , covering LSR velocities of 48-87 km s _1 
in G10.47+0.03 and 27-93 km s" 1 in SgrB2. The continuum 
was fitted using line-free channels (excluding 58-77, 44-80, and 
42-82 km s 1 in G10.47+0.03, SgrB2-N, and SgrB2-M, respec- 
tively) and subtracted from the line data. Imaging was done with 
almost uniform weighting for continuum maps (MIRIAD robust 
parameter -2, resulting in beams of 0.12"x0.06" forG10.47 and 
0.12" x 0.08" for SgrB2) and more natural weighting for line 
maps (robust 0.5, resulting in beams of 0.14" x 0.1" for G10.47 
and 0.15" x 0.11" for SgrB2). For total fluxes and lower noise 
levels, the continuum was also imaged in natural weighting. 

A pixel size of 0.02" and a cutoff for cleaning of twice the 
rms noise in the image was used. TableQjsummarizes beam sizes 
and noise levels. 



1 The National Radio Astronomy Observatory is a facility of the 
National Science Foundation operated under cooperative agreement by 
Associated Universities, Inc. 

2 http://bima.astro.umd.edu/miriad 



Table 3. Results of Gaussian fits to the emission lines shown in 

Fig.m 



field 


line width" 


velocity 


intensity 




(km s" 1 ) 


(km s" 1 ) 


(K) 


G10.47+0.03 south-Bl 


9.7 


69.4 


200 


G10.47+0.03 west-Bl 


10.8 


67.1 


170 


SgrB2-N east-K2 


14.9 


65.9 


100 


SgrB2-N west-K2 


22.9 


67.3 


100 


SgrB2-N south-K2 


11.9 


62.9 


100 


SgrB2-M north-F3 


6.8 


62.4 


100 



Notes. <n) Full Width at Half Maximum (FWHM) of the Gaussian 



For G10.47+0.03, the continuum clean components were 
used to self-calibrate the phases and apply the solutions to the 
line data. The absorption lines in this source were almost invis- 
ible before self-calibration. The procedure did not improve the 
SgrB2 data, presumably due to the good calibration with SgrA* 
only 0.7° away. All figures shown in this paper were created with 
the GILDAS software^. 

2.2. Effelsberg 

With the Effelsberg 100-m telescope of the Max-Planck-Institut 
fur Radioastronomie , we observed sev eral sources in HCN direct 
Mype lines in 2000 dThorwirthl200lh . 2002, and 2007. The lines 
were the 7=9 transition at 20.18 GHz, 7=10 at 24.66 GHz, 7=1 1 
at 29.58 GHz, and 7= 12 at 34.95 GHz. Beam sizes are 37,31, 26, 
and 22", respectively. Pointing was done on nearby quasars, and 
is accurate to about 10". The different polarizations and scans 
were averaged, and a baseline was subtracted by fitting a poly- 
nomial of order to 5 to line-free channels. The noise tube units 
measured by the telescope were converted to fluxes by compar- 
ison w ith a source of known flux, mostly NGC7027 dOtt et alj 
1994), which was observed once per day. To correct the defor- 
mation of the telescope at low elevations, the resulting flux was 
multiplied by an ad ditional factor of (1 - 0.01 x (30° - elv))~ l 
dGalli more et a l.l200ll) . which is 1.25 at an elevation of 10°, typ- 
ical for SgrB2. 

The observations were complemented by data from the 
millim eter-wave line survey of SgrB2 with the IRAM 30-m tele- 
scope dBelloche et alj|2008l) . 

3. Results 

For each source we show spectra (Fig. [TJ, continuum and 
integrated line maps, as well as position-velocity diagrams. 
Absorption lines are summarized in Tableland emission lines 
in Table El 

3 http://www.iram.fr/IRAMFR/GILDAS 
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Table 2. Parameters of absorption lines shown in Fig.Q] 



Hn region 


RA 


Dec. 


continuum" 


line width 6 


velocity 


line/continuum 


r 


HCN column density^ 




(J2000) 


(J2000) 


(K) 


(km s" 1 ) 


(km s" 1 ) 






(cm" 2 ) 


G10.47+0.03 Bl 


18 08 38.243 


-19 51 50.22 


2600 


18.9 


65.8 


0.72 


0.33 


5.5(19) 


G10.47+0.03 B2 


18 08 38.22 


-19 51 50.62 


1800 


18.5 


66.4 


0.776 


0.25 


4.1(19) 


G10.47+0.03 A 


18 08 38.188 


-19 51 49.74 


750 


10.3 


68 


0.67 


0.4 


3.6(19) 


SgrB2-N K2 


17 47 19.882 


-28 22 18.45 


1150 


10.0 


63.6 


0.16 


1.83 


1.6(20) 


SgrB2-N K3 


17 47 19.902 


-28 22 17.18 


1200 


7.3 


64.3 


0.652 


0.43 


2.8(19) 


SgrB2-MFlc 


17 47 20.126 


-28 23 03.9 


2650 


11.5 


60.4 


0.664 


0.41 


4.1(19) 


SgrB2-MFle 


17 47 20.131 


-28 23 04.02 


2300 


9.2 


60.9 


0.474 


0.75 


6.1(19) 


SgrB2-MFlf 


17 47 20.147 


-28 23 03.88 


3500 


15.3 


55.2 


0.714 


0.34 


4.4(19) 


SgrB2-M F3 


17 47 20.172 


-28 23 04.58 


5350 


6.7 


71.2 


0.696 


0.36 


2.1(19) 



Notes. The beam size is 0.14" x 0.1" for G10.47+0.03 and 0.15" x 0.11" for SgrB2. The rms on 1 km s~' channels is 120-150 K. Results of 
Gaussian fits and lower limits on the optical depth r and on the column density of hot HCN are given. The continuum flux density is the 
brightness temperature at the given locations in the line maps (which have a larger beam than the pure continuum maps). To convert from K to 
mJy/Beam, one has to divide by 53 for G10.47+0.03 and by 45 for SgrB2. {b) Full Width at Half Maximum (FWHM) of the Gaussian (c) The 
numbers in parentheses are powers of 10. 




18 > 0B m 38?30 38f25 38?20 3B?15 18 , 08'"38f30 38^5 3S?20 38?15 



Fig. 2. VLA maps of G10.47+0.03. The 7mm continuum emission is shown on the left, the integrated line map on the right (beams 
are depicted in the lower left). Contours denote 3 and 6 cr levels of the naturally weighted continuum (white, left) and of the inte- 
grated line emission convolved to 0.25" resolution (black, right). The white crosses give the positions toward which the absorption 
spectra were obtained and the white boxes the area over which the emission spectra are integrated, shown in Fig. [TJ The arrow 
denotes the direction along which the position-velocity plot of Fig. [3] was obtained. The map size is 2.5", about 0.13 pc. 



From Gaussian fits to the absorption lines, lower limits on 
the optical depth r and on the column density of hot HCN can 
be derived: The line-to-continuum ratio is <?~ T if line emission 
can be neglected, i.e. if the beam-averaged (excitation) tempera- 
ture of the absorbing gas is much lower than the beam-averaged 
background brightness temperature - else the line-to-continuum 
ratio is larger than e~ T . The background radiation is provided 
by free-free emission from ionized gas of around 10 4 K, which 
is at least an order of magnitude larger than the temperature of 
the absorbing gas. However, the observed continuum is lower 
than 10 4 K due to optically thi n free-f ree radiation in parts of the 
beam (see e.g. [Cesaroni et all (1201 Oh for spectral energy distri- 
butions of Hn regions with density gradient). Additionally, the 
hot molecular gas could emit over the whole beam. Therefore, 
the real optical depths can be much higher than derived from the 
line-to-continuum ratio, and the values given in Table|2]are only 
lower limits. 



From r and the observed line width Av (FWHM), the HCN 
column density can be derived, assuming Local Thermodynamic 
Equilibrium (LTE, see also Sect. 14. Il l at a temperature T: Using 



8ttv 2 



N u A u] <S>(e& -l) 



(1) 



with frequency v — 40.7669 GHz, Einstein A coefficient A u \ = 
3.75 x 10~ 8 s _1 and profile function at Gaussian line center <D = 

■J^^ x j; ~ ^r, the total column density N can be obtained 
from the column density in the upper state 



N u =N 



Q(T) 



(2) 



inserting the statistical weight g u = 2J + 1 = 27, the energy of 
the upper level E u = 1413 K, and the partition function Q(T) 
(Bar ber et al.l l2002). The temperature has two competing effects 
on these equations: While a higher temperature raises N u , it low- 
ers the population difference between upper and lower level (last 
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Fig. 1. The V2=l, 7=13 direct ^-type line of HCN observed to- 
ward several positions in G 10.47 +0.03 (left column), SgrB2-N 
(central column), and SgrB2-M (right column). The emission 
lines are averaged over the areas shown in Figs. |2] 2] and |6] 
while the absorption lines are observed with the beams given 
in Table Q] Overlaid in green are Gaussian fits, whose parame- 
ters are given in Tables [2] and [3] The vertical dashed lines mark 
the assumed systemic velocities. 



factor of Equation [TJ. For given r and Av, a minimum of N is 
reached around 500 K, while twice as many molecules are re- 
quired at around 300 and 1000 K, and 10 times as many at 200 
and 1800 K (100 times at 140 K). In Tabled we give this min- 
imum, which can be regarded as a lower limit on the column 
density of hot HCN. Only if the level population difference is 
higher than expected from LTE, the true HCN column density 
would be lower. 



3.1. G10.47+0.03 

The 7mm continuum (Fig. [2] left) shows the two hypercompact 
Hn regions Bl and B2 as well as the cometary-s haped ultracom- 
pact Hn region A (see also ICesaroni et alJ2010h . The total fluxes 
are 80 mJy for Bl, 48 mJy for B2, and 79 mJy for A. 

The integrated line map (Fig. [2] right) shows the absorp- 
tion and the distribution of line emission, which is concentrated 
around Bl and especially strong between Bl and B2. The ab- 
sorption lines toward Bl and B2 have a component at lower ve- 



T _ _ _ " - 1 



S 70 



IT-*" f^**^ 




G 
3 



-400 



Position (arcsec) 

Fig. 3. Position-velocity diagram of G10.47+0.03, along the ar- 
row shown in Fig. [2] The velocity of 68 km s _1 is marked as a 
dashed line, and the color scale is centered at K. 




■ 



200 



| 

a 
a 



PI, 



Position (arcsec) 

Fig. 5. Position-velocity diagram of SgrB2-N, along the arrow 
shown in Fig. [4] The velocity of 64 km s _1 is marked as a dashed 
line, and the color scale is centered at K. 



locities than the systemic velocity of around 68 km s 1 (derived 
from emission lines), indicating expansion motions (Figs. Q] and 
ED- 



3.2. SgrB2-N 

The c ontinuum map (Fig.|4]left) shows K2 and K3 dGaume et alJ 
1995) with fluxes of 70 and 140 mJy, respectively. Outside of the 
range shown in Fig. |4j we also detect K4, Kl (as a north-south 
elongated, 0.5" long Hn region of around 100 mJy), and a 10 
mJy point source at R.A. 17:47:19.905, Dec. -28:22:13.46 (4" 
north of K3). This latter source is associated with C2H5CN emis- 
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Fig. 4. VLA maps of SgrB2-N. The 7mm continuum emission is shown on the left, the integrated line map on the right (beams are 
depicted in the lower left). Contours denote 3 and 6 cr levels of the naturally weighted continuum (white, left) and of the integrated 
line emission convolved to 0.25" resolution (black, right). The white crosses give the positions toward which the absorption spectra 
were obtained and the white boxes the area over which the emission spectra are integrated, shown in Fig. [1] The arrow denotes the 
direction along which the position-velocity plot of Fig. was obtained. The map size is 3.5", about 0.13 pc. 



sionat 75 km s 1 and 3mm continuum detected bv lLiu & Snvderl 
(Il999h . We note th at emission from co mplex molecules (such as 
amino acetonitrile, iBelloche et al.l l2008) peaks at K2. 

The integrated line map (Fig.[4]right) shows wide-spread line 
emission around K2. Toward K2 itself, the line absorbs over 80% 
of the continuum, implying an enormous optical depth (Fig. Q] 
and Table|2|. The southern part seems to have a bit lower veloc- 
ities than the northern part (Figs. Q~]and |5). 

3.3. SgrB2-M 

The continuum map (Fig.[6]left) resembles the one obtained by 
Ide Pree et al.1 dl998l) . This source consists of a whole cluster of 
hypercompact Hn regions, including F3 with 1.3 Jy and Flc/e/f 
with about 0.5 Jy. 

Line absorption is concentrated toward the strongest sources 
Fl and F3 (Fig. [6] right). Toward Fl, the line is at 60 km s _I 
with a component at 50 km s _1 in Flf (Fig. [7] left). A northeast- 
southwest velocity gradient from 69 to 73 km s over 0.3" 
(2300 AU) is detected in F3 (Fig.Qright). In addition to the lines 
from Fig.Q] there is a possible absorption line toward F2 at 57 
km s" 1 with a width of 7 km s _1 and a line-to-continuum ratio of 
0.8. Although not visible in the integrated line map, there is also 
line emission present, as seen by averaging over a field outside 
the Hn regions. 

3.4. Single-dish data 

We investigated several sources with the Effelsberg 100-m tele- 
scope for HCN direct ^-type transitions, with no detections 
in G31.41+0.31, G34.26+0.15, W3(H 2 0), and W51d. Toward 
W51e, the 7=12 line was detected (with half the intensity of 
G10.47+0.03), but the lower-/ transitions are questionable (7=9 
at least 10 times lower than in G10.47+0.03). Clear detections 
in Galactic star-forming regions were only made in Orion-KL, 
G10.47+0.03, SgrB2-N, and -M. These latter three sources are 
shown in Fig. [8] with additional IRAM 30-m data for SgrB2. 



7=19 is contaminated by an Hy recombination line whose fre- 
quency corresponds to a ~20 km s higher velocity. 7=20 is 
blended with a transition of ethanol in SgrB2-N and possibly 
with Hey in -M. 

To compare our VLA data to the single-dish data, Fig. [9] 
shows the line integrated over the whole map of Figs. [2] |4] and|6] 
The VLA flux in G10.47+0.03 and SgrB2-N agrees well with the 
expectations from the Effelsberg data. In SgrB2-M, absorption 
dominates at low frequencies, but in the IRAM 30-m data weak 
emission can be seen. Emission also contributes to the VLA flux. 

A bit puzzling is the different widths of 7=9 and 10 in the 
Effelsberg data of SgrB2-M: It is 20 km s -1 in 7=9 and only 6.6 
km s _1 in 7=10 (Fig. [8}. As we found no blending at the 7=9 
frequency, this is probably due to different ratios of the 50, 60, 
and 70 km s _1 components: For 7=10, the 60 km s _1 component 
dominates. 

In SgrB2-N, also the higher-7 transitions with level energies 
up to 2100 K were clearly detected. This allows to construct a 
rotation diagram (Fig.fTOli. yielding a temperature of 485 + 50 K 
and an HCN column density in a 1" source of (5.8 + 2) x 10 19 
cirT 2 . Note that taking absorption and optical depth into account 
would lead to a higher column density and a lower temperature, 
since both effects weaken especially the lower-7 lines. 



4. Modeling 

In an attempt to constrain the spatial structure of the hot molecu- 
lar gas, we constructed radiative-transfer models that reproduce 
the observations. We employed a trial-and-error technique and 
compared the model to the data by eye. This method yields (at 
best) a model that is consistent with the data, but it cannot pro- 
vide errors of the parameters, and different models could fit as 
well. The three-dimensional radiative-transfer code RADMC- 
3E0 was used to compute the radiation that a model emits. It al- 
lows adaptive mesh refinement, following user-defined criteria. 



4 http://www.ita.uni-heidelberg.de/~dunemond/software/radmc-3d 
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Fig. 6. VLA maps of SgrB2-M. The 7mm continuum emission is shown on the left, the integrated line map on the right (beams are 
depicted in the lower left). Contours denote 3 and 6 cr levels of the naturally weighted continuum (white, left) and of the integrated 
line emission convolved to 0.25" resolution (black, right). The white crosses give the positions toward which the absorption spectra 
were obtained and the white boxes the area over which the emission spectra are integrated, shown in Fig. [1] The arrow denotes the 
direction along which the position-velocity plot of Fig. |7]was obtained. The map size is 4", about 0. 15 pc. 




Position (arcsec) Position (arcsec) 

Fig. 7. Position-velocity diagrams of SgrB2-M, along the arrows shown in Fig. [6] The cut through Fl is shown on the left, F3 on the 
right. The velocity of 60 km s _1 is marked as a dashed line, and the color scale is dominated by absorption. 



Dust properties and the dust density distribution as well as loca- 
tion, surface temperature and radius of stars are given as input. 
The program then computes the dust temperature by tracing pho- 
ton packets which are randomly emitted by the stars. Additional 
inputs are density and temperature of ionized gas, the velocity 
field, fundamental molecular data, and the molecular abundance. 

4.1. Assumptions 

We use the dust opacity from Ibssenkop f & Henningl (Il994l) 
without grain mantles or coagulation. This is reasonable assum- 
ing that ice mantles around dust grains have evaporated com- 
pletely and recently, so that the grains have had no time to coag- 



ulate again. Since the VLA continuum is dominated by free-free 
emission, the data are not sensitive to the dust emission and so 
we cannot constrain the dust optical depth, which is a combi- 
nation of dust opacity and density, through the dust emission. 
However, the heating is strongly affected by the dust optical 
depth. 

The temperature of the ionized gas is assumed to be 10 4 K, 
which is the order of magnitude expected from cooling by trace 
species. A different temperature would require a different elec- 
tron density to account for the observed fluxes. 

The HCN molecular data are from iThorwirth et al.l 
(2003a) and ta ken from the Cologne Database for Molecular 
Spectroscopy dMuller et al.ll200ll l2005h . We use an HCN abun- 
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G1 0.47+0.03 SgrB2-N SgrB2-M 




20 40 60 80 100 20 40 60 80 100 

Velocity (km/s) Velocity (km/s) 

Fig. 8. Single-dish spectra of G 10.47 +0.03 (left), SgrB2-N (cen- 
tral), and SgrB2-M (right). 7=9 (20 GHz) to 12 (35 GHz) were 
observed with the Effelsberg 100-m telescope, J =19 (85 GHz) 
to 22 (113 GHz) with the IRAM 30-m telescope. Beam sizes 
vary between 20 and 37", so the whole flux is received by the 
telescope. 



G1 0.47+0.03 SgrB2-N SgrB2-M 
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Velocity (km/s) Velocity (km/s) Velocity (km/s) 



Fig. 9. VLA line flux of G10.47+0.03 (left), SgrB2-N (central), 
and SgrB2-M (right), integrated over the whole maps shown in 
Figs. GO ED andil 



dance of 10~ 5 relative to H2 for our models. Since we mainly 
constrain the density of HCN, a different abundance would re- 
quire different dust densities (with the gas/dust mass ratio of 
100). Dust continuum data (from the SMA) indicate that the re- 
sulting dust densities are on the right order of magnitude, which 
justifies this abundance assumption. 

The line radiative transfer assumes LTE (full non-LTE radia- 
tive transfer is also planned for RADMC-3D). LTE is a reason- 
able assumption for the observed line as HCN thermalizes to the 
ambient dust temperature: It is vibrationally excited by 14jt/m 




1000 1500 2000 

E upp „ / k (K) 



Fig. 10. Rotation diagram for SgrB2-N, based on the single-dish 
fluxes. The integrated line flux is source-averaged for a source 
size of 1". The first four data points are from Effelsberg, the 
last four from IRAM 30-m; the VLA flux fits very well to the 
expectation. 



radiation emitted by warm dust, which is very optically thick at 
this wavelength due to the high densities (an optical depth of 1 
is reached after about 25 AU in a density of 10 8 cm -3 ). In this 
warm and dense environment, also the levels in the ground vi- 
brational state are thermalized by infrared pumping as well as by 
collisions and radiative excitation. Deviations from LTE occur at 
low temperatures, where the vibrational levels are not populated 
anyway. 

The models have no macroscopic velocity field. We assume 
a constant intrinsic line width (microturbulence) throughout the 
source, which reflects internal motions. The resulting spectra are 
shifted by the source velocity (68 km s _1 for G10.47+0.03, 64 
km s" 1 for SgrB2-N, 61 km s" 1 for SgrB2-M Fie, and 71 km s" 1 
for SgrB2-M F3). 

4.2. Modeling results 

To compare the model to the observational data, the synthetic 
map produced by RADMC-3D, supplemented by distance and 
coordinate information, is Fourier- transformed, folded with the 
uv coverage of the observations, and imaged again. The first 
step of our fitting procedure is to reproduce the observed contin- 
uum by creating H n regio ns whose sizes and densi ties are based 
on ICesaroni et al.l (1201 Oh for G 10.47 +0.03 and Ide Pree et all 
dl998h for SgrB2-M. For SgrB2-N, we u se similar models as 
for G10.47+0.03 since the 1.3cm data of lGaume et"al] d 1995b 
are not sufficient in terms of angular resolution and wavelength. 
Sizes and densities of the Hn regions are slightly modified to fit 
the observed continuum. The hypercompact Hn regions Bl and 
B2 in G10.47+0.03 and K2 in SgrB2-N have thus density gradi- 
ents (with r -15 and outer radii of 1590 AU in Bl, 950 AU in B2 
and 1000 AU in K2), while all other Hn regions have constant 
densities. 

ICesaroni et all d2010l) and Ide Pree et all (1 19981) also give 
spectral types of the excit i ng star s, which are converted to lumi- 
nosity following Panagiaj (Il973h . These stars heat the dust and 
are, for simplicity, placed in the plane of the sky, except for two 
Hn regions in SgrB2-M. The dust density distribution is adapted 
to fit the line observations. FigureQT]shows the column densities 
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in x, y, and z directions of the three models that are described in 
the following. 

G1 0.47+0.03 The model consists of a clump centered at Bl 
whose density follows a Gaussian with 7 x 10 7 H2 ctrT 3 at the 
half-power radius of 7000 AU. This dust is heated by Bl with 
10 5 L , B2 with 8.3 x 10 4 L , and A with 4.6 x 10 4 L . The 
intrinsic line width (FWHM) is 8.3 km s -1 . The model output is 
compared to the observations in Fig. [12] 

SgrB2-N In the model for SgrB2-N, the dust is heated by K2 
with 10 5 L and K3 with 8.3 x 10 4 L . The dust density follows 
a Gaussian centered at K2 with 7 x 10 7 H2 cirT 3 at the half -power 
radius of 10 4 AU. An additional core (Gaussian with 5 x 10 8 H2 
cm -3 at the half-power radius of 1000 AU) 3000 AU in front of 
K2 provides the strong absorption. The intrinsic line width is 10 
km s . Figure [T3l shows the comparison to the observations. 

SgrB2-M In this model, the dust density follows two radial 
power laws with index 1.5, centered at Flf and F3. The first one 
starts with 5.5 x 10 8 H2 crrT 3 at a radius of 660 AU, the latter one 
with 5 x 10 8 H 2 cirT 3 at 750 AU. Inside this latter radius, there 
is an Hn region that contains a star of 8 x 10 4 L . Additionally, 
a larger Hn region is located 5000 AU closer to us with a star of 
2.5 x 10 5 L Q - F3 is thus split into these two different Hn regions. 
Fie with 2.5 x 10 4 L is located 2000 AU farther away from 
us than the other Hn regions to facilitate absorption (Fig. ITTb . 
Further heating is provided by Flc, Flf, F2, F4 (each 5.4 x 10 4 
L©), Fla (3.8 x 10 4 L ), and F3a (2.5 x 10 4 L ). The intrinsic line 
width is 6.7 km s _1 . A comparison to the observations is shown 
in Fig. [Ml 

These models are as simple as possible to approach a good fit 
to the observations, but are not unique. The model maps shown 
in Figs. [12] to [14] would appear more similar to the data if cor- 
responding noise was added. This noise also makes it difficult 
to constrain more complicated models, since it is not clear how 
much of the clumpy substructure is real (which is on the order 
of only 2 or 3 times the rms noise) 

5. Discussion 

In this section, we discuss implications from the observational 
results and the modeling efforts. 

5. 1 . Optical depth 

As explained in Sect. [3] lower limits on optical depth and 
HCN column density can be derived from the absorption lines 
(Table [2]). Outstanding is K2 in SgrB2-N, where the line is op- 
tically thick. The column densities of hot HCN toward the Hn 
regions where it is detected are at least a few times 10 19 cm -2 , 
which translates into H2 column densities of 10 24 to 10 25 cirT 2 
even for a high HCN fractional abundance of 10~ 5 . This rep- 
resents only the molecular gas at temperatures above roughly 
300 K. In the models, the mass above 300 K is 350 Mo in 
G10.47+0.03, 500 M in SgrB2-N, and 200 M in SgrB2-M. 

Such large quantities of hot molecular gas are necessary to 
explain the observations, and challenge theories of massive star 
formation. Current hydrodyn amical simu lations produce far too 
little hot molecular gas (e.g. lPeters et al.ll2010L with a total mass 
of 1000 Mo), both absolute and relative to the total mass, since 
the highly inhomogeneous structure (filamentary disks) does not 



lead to diffusion of the heating radiation. This is probably due to 
the lower initial mass and hence lower column densities and dust 
optical depths of these models, and it is to be hoped that in the 
future simulations of cores as massive as the ones presented in 
the current paper will reproduce our results. Also a high abun- 
dance of HCN (on the order of 10~ 5 ) in the dense, warm gas 
is needed, since the dust densities of our models are consistent 
with dust continuum data from the SMA and a lower abundance 
would mean an even higher mass of the hot gas. We expect fu- 
ture chemical models to compute such a high HCN abundance, 
probably through high-temperature gas-phase reactions. 

Although HCN thermalizes to the ambient dust temperature, 
the population difference between the two {-type levels is very 
small, making them sensitive to the exact excitation conditions. 
The line optical depth is proportional to the population differ- 
ence, 

1411 K 1413 K 

t oc N\ - N u oc e T > - e ^ , (3) 

where T\ and T u are the excitation temperatures for the lower 
and upper state, respectively. A difference of 0.1% between T\ 
and T u causes t to change by more than 70%; inversion occurs 
if T u is 0.14% larger than T\. The two levels are independently 
coupled to the ground vibrational state through the 14/vm tran- 
sitions, which determine the excitation temperatures. If these 
are different, the analysis would be far more difficult, as one 
would have to know the exact 14//m radiation field, where lines 
could slightly affect the optical depths. However, there is no such 
mechanism known, and an extreme sub-thermal population dif- 
ference, which would greatly reduce the required absorption col- 
umn densities, is highly unlikely because of the agreement be- 
tween absorption and emission, between 7=13 and other direct 
£-type lines (Fig. [Tol l, and between the derived column densi- 
ties and dust emission data. Therefore, we consider LTE a good 
assumption for this transition. 

If we assume a similar optical depth for the emission lines, 
say 0.3, and a temperature of 400 K, then the line intensity 
should be 100 K, which is just the noise level in 1 km s chan- 
nels and 0.13" beams. Owing to the strength of the free-free ra- 
diation, we are thus more sensitive to absorption than emission 
lines. Averaging over a larger region (Fig. [T]and Table |3} lowers 
both noise and peak intensities. The latter reach up to 300 K in 
G10.47+0.03, suggesting that optical depth and temperature are 
higher than the above mentioned assumptions. Surprising is the 
detection of the line towar d G10.47+0.03 A, sin ce no NH 3 (4,4) 
satellites were detected by Cesaron i et al.l (l2010l) . 

In the model for SgrB2-M, we have placed the Hn regions F3 
and Fie in front of and behind the bulk of the molecular gas, re- 
spectively. This prevents absorption toward most of F3 and leads 
to stronger absorption toward Fie. The exact offsets in z are not 
well constrained, however. To obtain the very strong absorption 
toward K2 in SgrB2-N, we put a very dense core just in front of 
the Hn region. This is of course not very satisfying as it is un- 
likely that such a core is exactly along the line-of-sight. It was 
however not possible to reproduce both emission and absorption 
with a more symmetric distribution. 

5.2. Heating 

Heating up large masses of molecular gas requires deeply em- 
bedded massive (proto)stars with high lumi nosities. Their radi- 
ation is either originally in the infrared (Hoso kawa & Omukail 
2009), not producing Hn regions, or is quickly processed to the 
infrared by dust absorption of the UV radiation. Due to high col- 
umn densities in all directions, the dust is optically thick even 



8 





G10.47+0.03 



SgrB2-N 



-1. 10* 




-1. 10* 



-1. Ifi 4 



y (au) 




y (au) 































































o 
















- - 











































x (MJ) 




-1. 10* Ox (AU)10 4 



-1. Ifi* Ox (MtflO* 



| i i i i | i i i i | 



\ I I I I I I I I I I I I I \ 

-1. 10* Ox (M010 4 



~— 4 id" 



■ 2 10* 



Fig. 11. Column density maps of the models for the three sources. The lower panels show the column density (in H2 cm -2 ) in 
z-direction, the central panels in y-direction, and the upper panels in x-direction (the arrow points to the observer, who is at positive 
z). Two Hn regions are marked for each source. The model radiation is compared to the observations in Figs. H2ltofT4l 



in the infrared and this radiation cannot escape, but diffuses out- 
wards by multiple absorption/emission events until the dust is 
optically thin to its own radiation. This diffusion (or radiative 
trapping) leads to much higher temperatures in the inner part; 
e.g. at 2000 AU from the HCH11 region Bl, diffusion raises the 
temperature from 150 K in a low-density model, where diffu- 
sion is not effective, to 450 K in the high-density model for 
G10.47+0.03. This mechanism therefore seems indispensable to 
account for large masses of hot gas. Additionally, by absorbing 
stellar UV radiation the dust protects the molecules from dissoci- 
ation. In general, heating is mo re efficient if the heati ng source is 
inside a density condensation dKaufman et alj[l998f) than if it is 
external, as could be the cas e for G34.26+0.15 (IWatt& Mundv 
ll999UMookeriea et al.ll2007h . 

In our models, the diffusion mechanism is accomplished by 
dense spherical clumps around heating sources. It is however 
not clear how many heating sources there are: The stars in the 
Hn regions are sufficient to reproduce the observed line strength, 
as the presented models demonstrate. But a good fit can also be 
achieved when lower-luminosity stars are added to the models 
according to the Initial Luminosity Function, without changing 
the stars in the Hn regions. So on the basis of these data it seems 
not possible to distinguish between the scenarios of heating only 
by the stars ionizing the Hn regions and heating by Hn regions 



plus embedded sources which do not emit in the radio regime. 
We may speculate that there are much more stars present than 
detected from the free-free radiation, since the observed sources 
are clusters in formation. Also, there could be more extended 
(faint) line emission (and the absorption coming from a larger 
column) than in the models, which would point to additional 
heating sources. Furthermore, the simple models predict strong 
self-absorption in the (optically thick) submm rotational transi- 
tions within the V2=l state of HCN and H 13 CN, which is how- 
ever not observed with the SMA, presumably due to dumpiness 
(Rolffs et al. 201 1, in prep.). 



While not excluding the possibility of additional heating 
sources, we could show that heating by the stars in the hyper- 
compact Hn regions is sufficient to reproduce the observations. 
Hence, our modeling results suggest that it is unlikely that the 
main heating is due to embedded objects without ionized gas. 
Given that the sources are dense condensations deeply embed- 
ded in giant molecular clouds and have high column densities, 
diffusion of radiation, which is necessary to increase the heat- 
ing, is reasonable even for more inhomogeneous and asymmetric 
structures than modeled. 
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Fig. 12. Model for G10.47+0.03 (right) compared to the data (left). The model spectra are also overlaid in red. The white contours 
denote the continuum in steps of 1000 K. 



5.3. Velocity field 

Many systematic motions are hidden in the line width of around 
10 km s _1 , which contains both turbulent and macroscopic mo- 
tions. These components (e.g. the presence of rotation, infall, 
or outflows) can not be disentangled. For all emission lines, the 
signal-to-noise ratio is too poor to gain accurate information on 
the velocity structure. 



G1 0.47+0.03 The most striking feature in this source is an ad- 
ditional blue-shifted absorption component in both Bl and B2 at 
a velocity of around 60 km s _1 , compared to the source velocity 
of 68 km s" 1 (Fig Q]i. The NH 3 (4,4) absorption (ICesaroni et al.1 
rT998ll20Toh is at 53 km s 1 toward Bl and at 61 km s 1 toward 
B2, with much lower optical depth than the emission compo- 
nent. An explanation could be that the colder gas, probed by 
NH3(4,4) with its level energy of 200 K, is mainly in an outflow 
cone, while the hot gas probed by HCN,V2=1,./=13 is mainly at 
the systemic velocity close to the HCHu regions and only partly 
expanding. With the SMA, we see broad blue-shifted absorp- 
tion in several lines toward the dust continuum peak between 
Bl and B2 (Rolffs et al., in prep.). This indicates expansion mo- 
tions involving many different molecular excitation c onditions, 
while at the same time large-scale infall is present (Rolf fs et al.l 
1201 lb . Given the high excitation and the different velocities, an 
unrelated foreground velocity component is unlikely. The region 
west of Bl appears to move at a velocity of 67 km s _1 , while 
south of B 1 the velocity is 69 km s~ 1 . 



SgrB2-N While the strong absorption line toward K2 has a line 
width of 10 km s _1 at the source velocity of 64 km s _1 , the emis- 



sion lines north of K2 have a higher line width and higher veloc- 
ity, possibly indicating an additional component at ~70 km s . 

SgrB2-M Three velocity components can be seen from the ab- 
sorption lines (Fig. [7]i, one at around 50 km s in Flf, one at 60 
km s _1 in Flc, Fie, and Flf, and one at around 70 km s _I in F3. 
Relative to the ambient velocity of 60 or 64 km/s, which itself 
is not well determined due to different components, the Flf line 
shows blue-shifted absorption, indicating an expansion compo- 
nent, and the F3 line is red-shifted. This could mean infall or a 
motion of the whole F3 core. The F3 line shows a northeast- 
southwest velocity gradient from 69 to 73 km s over 0.3" 
(2300 AU). If that is interpreted as Keplerian rotation, and if it 
is assumed that the whole orbit is seen and is edge-on, the lower 
limit on the central mass would be 5 M . The large-scale veloc- 
ity fie ld shows expansio n in the inner parts and infall in the outer 
parts (iRolffs et alj|2"oloh . 



6. Conclusions 

With the VLA, we have obtained maps of a high-excitation 
(1400 K) HCN line at high angular resolution (0.1") toward 
three massive star-forming regions. We see both absorption to- 
ward small and dense Hn regions and emission. The infered op- 
tical depths require on the order of 10 24 H2 cirT 2 and 100 solar 
masses of hot (>300 K) molecular gas within a region of less 
than 0. 1 pc. We note that the column densities of hot HCN de- 
rived from the absorption lines are only lower limits, and can be 
substantially higher especially due to beam dilution. Heating is 
provided by the stars powering hypercompact Hn regions, pos- 
sibly complemented by heating sources invisible in the radio 
domain. Diffusion of radiation in an optically thick dusty en- 
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VIA HCN.v 2 -1 J=13 




Fig. 13. Model for SgrB2-N (right) compared to the data (left), 
denote the continuum in steps of 1000 K. 

vironment is likely to increase the temperatures. The line ve- 
locities reveal expansion motions in G10.47+0.03 and compo- 
nents at 50, 60, and 70 km s -1 in SgrB2-M, the latter with a ve- 
locity gradient. Our three-dimensional radiative transfer model- 
ing with RADMC-3D demonstrate the great power and potential 
of our approach to constrain the source structure of high-mass 
star-forming regions. HCN direct ^-type lines are a unique tool, 
whose future observation will also be possible with ALMA and 
the SKA. 
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VIA HCN.v 2 =1 J=13 model for SgrB2-M 




Fig. 14. Model for SgrB2-M (right) compared to the data (left). The model spectra are also overlaid in red. The white contours 
denote the continuum in steps of 1000 K. 
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